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Abstract 

A generalization of the core/halo model to the Bose-Einstein correlation function 
of any number of particles is presented. In particular, a simple prediction is obtained 
CN ' for the intercept parameter of the n-particle Bose-Einstein correlation functions for 

> ; 

fS| arbitrary large values of n. 

in 
o 

Introduction. Multi-particle correlations are studied extensively in high energy colli- 
sions, with the aim of learning more about possible phase transitions or self-similar fluctu- 
ations. Especially, the short-range part of the two-particle correlation function, is thought 
to carry information about the space-time structure based on a quantum-statistical effect 
discovered in refs. [1], Q. One of the most important driving force behind correlation 
studies in high energy physics is the possibility to measure space-time dimensions on the 
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10~ 15 m and 10~ 23 sec scale. See refs. || |j] for recent reviews on high energy heavy ion 
physics, for summaries on correlation studies in this field see refs. [|, Q [7j, 

Multi-particle correlations are becoming experimentally well determinable in present 
and planned heavy ion reactions at at Brookhaven AGS, at CERN SPS, at the Relativistic 
Heavy Ion Collider (RHIC) and at CERN LHC both due to the very large number of the 
produced particles and due to the dramatic increase of the data quality accessible already 
for multidimensional correlation studies, ref. ||. 

There has been a lot of development recently in the theory of multi-particle Bose- 
Einstein correlations, form partial coherence to multi-particle wave-packet models and 
event generators with multi-particle symmetrization, refs. JCD, O, |T2|, |TB|, [TJ, [T5|, [R| [T7|] . 



The purpose of the present Letter is to study the structure of multi-particle Bose-Einstein 
correlations in one special case, when the emission function can be well separated to a 
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core and a halo part. Such a scenario is referred to as the core/halo model The 



formal study of the core/halo type of models first started with numerical simulations, [|T^, 

m m m, p 



f23| , P4[ . The essential ideas necessary to formulate the core/halo model were 
also published in ref. [[ZIJ, see refs. |26| , |2T| for further details. 

In refs. |f2"B| , E^ , |30|j an Edgeworth expansion method was proposed to characterize the 
non-Gaussian features of the Bose-Einstein correlation functions as caused by interference 
of decay product pions from long-lived resonances. However, the data indicated surpris- 
ingly regular, simple Gaussian correlation functions when the analysis was performed in 
two or three dimensions pT|, |32|. This resulted |33[ in the re-formulation and simplification 
of the theoretical aspects of the resonance-decay effects on the Bose-Einstein correlation 
functions in a particularly transparent form, with an emphasis on the interpretation and 
possible utilization of the momentum dependent intercept parameter of the two particle 



Bose-Einstein correlation function in ref. |]18[ . We present herewith the generalization of 
this core/halo model to the n-particle case, utilizing the Wigner-function formalism. 

In this formalism, the one-boson emission is characterized by the emission function, 
S(x,p). Here x = (t,r) and p = (E, p ) denote four- vectors in space-time and in momen- 
tum space, for particles with mass m = y/ E' 2 — p 2 . An auxiliary quantity is 



S{Ak,K) 



d A x S(x, K) exp(iAA; • x) 



where Ak = pi — P2, K = (pi +p2)/2 and Ak ■ x stands for the inner product of the 
four-vectors. The invariant momenum distribution reads as 

dn 



E- 



dp 



N x (p) = S(Ak = 0,K = p). 



(2) 



In this Letter we shall utilize the 'hydrodynamical' normalization of the Wigner-functions, 



dp 
~E 



d A xS(x,p) 



dp 

~E 



Ni{p) 



(n), 



(3) 



where (n) is the mean multiplicity. The two-particle BECF-s are prescribed as 

l S(Ak,K)\ 2 \S(Ak,K)\ 2 



C2(P2,P 2 ) 



1 + — 



1 + 



S{0,K) 



|2 ' 



(4) 



s(o, Pl )s(o,p 2 r 

as was presented e.g. in ref. |2(], [7|. The precision of the last approximation was estimated 
to be of the order of 5 % as discussed in Ref. ||34|| . In this Letter final state interactions 
are neglected and a completely chaotic particle emission is assumed. 

One can show that the higher order Bose-Einstein correlation functions, 

N n (Pl, P2, Pn 



Cn(Pl,P2, 



iV 1 (p 1 )iV 1 (p 2 )...iV 1 (p n ) 



(5) 
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are given in terms of the Fourier-transformed Wigner-functions as 

n 

zn^ 

i i=i 

n 



C n (pi,p 2 , ...,p„) = g " n 1 = Y.f[-^7-T- = £II S ( i ' <7 *)) ( 6 ) 



a n i=l 



i=l 



where the summation is over all the a n permutations of n indexes, 

S (i c 

S(i,ai) = S(K ij<Ti ,Aki iai ) 7 and s(i,ai) = ~ \ - , (7) 

S{i,i) 

Ki,ai = P% + 2 Pai , and Ak i>ai = Pi - p ai . (8) 

Note the distinction between a n , which stands for the set of the permutations of the 
numbers l,2,...,n and the crj (subscript i), which indicates the number replacing the 
index i in a given permutation from the set a n . Note also that s(i, &{) = s*(ffh z ') 

§*(<Ji,i), although the relationship S(i,<Ti) = S*(ai,i) is satisfied. 
In the core/halo model, the following assumptions are made : 

Assumption 0: The emission function does not have a no-scale, power-law like struc- 



ture. This possibility was discussed and related to intermittency in ref. [37 



Assumption 1: The bosons are emitted either from a central part or from the surround- 
ing halo. Their emission functions are indicated by S c (x,p) and Sh{x,p), respectively. 
According to this assumption, the complete emission function can be written as 

S(x,p) = S c (x,p) + S h (x,p), (9) 

using the hydrodynamic normalization of the Wigner functions. 

Assumption 2: We assume that the emission function which characterizes the halo 
changes on a scale Rh which is larger than R max ~ h~/Q m in, the maximum length-scale 
resolvable [Tj| by the intensity interferometry microscope. However, the smaller central 



part of size R c is assumed to be resolvable, Rh > Rmax > R c - This inequality is assumed 
to be satisfied by all characteristic scales in the halo and in the central part, e.g. in 
case the side, out or longitudinal components [3B, ^T|] of the correlation function are not 
identical. 

Assumption 3: The momentum-dependent core fraction f c (i) = N c (pi) / N\(pi) varies 
slowly on the relative momentum scale given by the correlator of the core s c (l, 2)s c (2, 1). 

According to this smoothness assumption, which was only implicitly made in ref. fT8| , 
the relative momentum dependence of the core fraction can be neglected on the relative 
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momentum scales on which the correlator of the core is changing. If a typical core radius 
parameter in a given direction is indicated by R c then the slow variation of f c (i) implies 

that / c (Pi)/c(pj) ^ fc( K i,j) if I P» ~ Pi I < h / R c- 

The emission function of the center and that of the halo are normalized here as 

J d 4 x-^-S c (x,p) = (n) c , and J d 4 x-^-S h (x,p) = (n) h , (10) 

where the subscripts c, h index the contributions by the central and the halo parts, re- 
spectively. In ref |18|], the above Assumptions were formulated with the help of the 
Wigher-functions normalized to 1, and a core-fraction f c = (n) c /(n) was used, while in 
this paper the momentum- dependent core fraction / c (p) is utilized. One finds that 

Nx(p) = N c (p)+N h (p), and (n) = (n) c + (n) h . (11) 

Note, that in principle the core as well as the halo part of the emission function could be 
decomposed into more detailed contributions. In case of pions and NA44 acceptance , 
one can separate the contribution of various long-lived resonances as 

S H (x,p) = £ S^ lo (x,p) and N h (p) = £ N^ (p). (12) 

r=u},rj,ri l r=uj,n,r]' ,Kg 

For our considerations, this separation is indifferent. According to our assumptions, the 
Fourier - transformed Wigner-functions characterizing the (full) halo part are narrower 
than the two-particle momentum resolution Q m i n , for which we typically may take the 
value of cca. 10 MeV. It is important to keep in mind that the halo part of the emission 
function is defined with respect to Qmin- For example, if Q m in = 10 — 15 MeV, the decay 



products of the uo resonances can be taken as parts of the halo |3q] , however, should 
the experimental resolution and the error bars on the measurable correlation funtion 
decrease significantly below h/T u = 8 MeV, the decay products of the u resonances 



would contribute to the core. See refs. |38], [1J| for a more complete discussion on this 
specific topic. 

The measured Sh{ii<Ji) is thus vanishing if i ^ o~i at the given measured relative 
momentum Q^ ai > Qmin- Note that Qi jCH can be one, two- or three-dimensional quantity, 
e.g. any reasonable measure of Aki )(Ti . However, it is important to observe that Sh{i,i) 
contributes to the spectrum, and is not affected by the two-particle momentum resolution: 

Sh(i,cn) = 6i 7(Ti S(i,i), (13) 
s(i,<Ti) = S i)a . + (l-S i)(Ti )f c (i)s c (i,cri), (14) 
fc(i) = /c(Pi) = #c(Pi)/iVi(Pi)- (15) 
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Thus the Bose-Einstein correlation function C„(l, 2, n) = C(pi, P2, p n ) reads as 

n 

C n (l,2,...n) = Ell^i + C 1 -^)/^)^^)] ( 16 ) 



Let us introduce p n to stand for those permutations of (l,...,n) which are mixing all 
the numbers from 1 to n and let us indicate by pi the value which is replaced by % in a 
given permutation belonging to the set of permutations p n . (Superscript indexes a set of 
permutations, subscript stands for a given value). Then we have 7^ % for all values of 
i — 1, ...,n, while Oi = i type of replacements are allowed. 
Thus the general expression for C n (l, n) reads as 



n n 1 



C n {l,...,n) = 1 + E E EII/c(W*.i»)- (17) 

j=2i 1 ,...,i j =l p n k=l 

Here 53' indicates that the summation should be taken over those set of values of the 
indices which do not contain any value more than once. 

Let us explicitly write out the above expression until the I = 4 terms: 

n 

C n (l,...,n) = 1+ E fc(i)fc(j)[sc{i,j)s c {j,i)} + 
i^j=i 

n 

+ E fc(i)fc(j)fc(k) [s c (i,j)s c (j,k)s c (k,i) + s c (i,k)s c (k,j)s c (j,i)] + 

ii=j¥=k=\ 

+ e mfc(j)f c (k)fs) x 

x [s c (i, j)s c (j, k)s c (k, l)s c (l, i) + s c (i,j)s c (j, l)s c (l, k)s c (k, i)+ 

+ s c (i,k)s c (k,j)s c (j,l)s c (l,i) + s c (i, k)s c (k,l)s c (l, j)s c (j,i) + 

+ s c (i,l)s c (l,j)s c (j, k)s c (k,i) + s c (i,l)s c (l,k)s c (k,j)s c (j,i) + 

+ s c {i,j)s c (j,i)s c (k,l)s c (l,k) + s c (i,k)s c (k,i)s c (j,l)s c (l,j) + 

+ s c (i,l)s c (l,i)s c (j,k)s c (k,j)} + (18) 

There are 44 (and 265) terms in the next part which fully mixes 5 (and 6) different indexes 
for any given fixed value of i 7^ j 7^ k 7^ / 7^ n), respectively. 



For the two-particle case, we recover the earlier result given in Ref. ||18|| . Let us proceed 
very carefully at this point, so that the role of Assumption 3 could be made transparent. 
If only two particles are present, the above equation reduces to 

C 2 (l,2) = l + / c (l)/ c (2)S c (l,2)s c (2,l). (19) 
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If Assumption 3 is also satisfied by some experimental data set, then the resulting formula 
becomes particularly simple: 



C(Ak 12 ,K 12 ) = l + K(K 12 ) \ S / Ak ^ 1 ~l + K(K 12 ) \f£ k ' K) J ,(20) 

S c (0 lPl )S c {0,p 2 ) \S C (0, K l2 )\ 2 

(21) 

where the effective intercept parameter A*(iT 12 ) is given as \*(K 12 ) = [N C (K 12 ) / N^K^)} 2 . 
As emphasized in Ref. |18|], this effective intercept parameter ( ^ exact intercept parameter 
at Q = MeV) shall in general depend on the mean momentum of the observed boson pair, 
which within the errors of Qmin coincides with any of the on-shell four-momentum p% or 
p 2 . Thus one obtains the core/halo interpretation of the two-particle correlation function: 
the intercept parameter \*(K) measures the momentum dependent core fraction and the 
relative momentum dependent part of the two-particle correlation function is determined 
by the core. Note, that in high energy heavy ion collisions the momentum dependence of 
the A*(iQ parameter is very weak, actually, within the errors \*(K) is constant for the 



NA44 data analyzed in ref. [FJ|. However, the validity of this Assumption 3 has to be 
checked experimentally for each data set, by determining the momentum dependence of 
the X(K) parameter of the two-particle correlation function. 

Within the core/halo picture, the n-particle correlation function has a simple form if 
all the n momenta are approximately equal, i.e. | A/cj )CTj |< Q m i n for all i ^ o~i, and this 
situation is denoted by pi ~ p 2 — ... — p n — P. One obtains that 

c n ( Pl ~ P2 ~ ... ~ Pn ~ p) = i +E/ C (p) j n A <*i 

j= i 

where <x, stands for the number of fully mixing permutations of j indexes, i.e. the number 
of permutations in pP . The first few values of aj are given as 

«i = 0, a 2 = 1, a 3 = 2, a 4 = 9, a 5 = 44, a e = 265. (22) 

These values can be obtained from a recurrence relation, as follows. Let us indicate 
the number of permutations that completely mix exactly j non-identical elements by 
ay. There are exactly ( J ^ different ways to choose j different elements from among 
n different elements. Since all the n\ permutations can be written as a sum over the 
fully mixing permutations, the counting rule yields n\ = 1 + Y%=i ( j \ a j? which can be 
rewritten as a recurrence relation for ctf 

a n = n\-l-J2 ("W (23) 
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We have the following explicit expressions for the first few intercept parameters: 



A*, 2 


;p) 


= fc(P) 2 






(24) 


A*, 3 


;p) 


= 3/ c (P) 2 


+ 2/ c (P) 3 




(25) 


A #) 4 


;p) 


= 6/ c (P) 2 


+ 8/ c (P) 3 - 


f 9/ c (P) 4 


(26) 


A*, 5 


;p) 


= 10/ C (P) 


2 + 2o/ c (p; 


3 + 45/ c (P) 4 + 44/ c (P) 5 


(27) 


A*,6 


;p) 


= 15/ C (P) 


2 + 4o/ c (p; 


3 + 135/ C (P) 4 + 264/ c (P) 5 + 265/ c (P) 6 


(28) 



In general, the intercept parameter of the n-particle correlation function reads as 

n , v 

A*,n(P) = E/eCPyf^)^ (29) 

where a.,- is defined by the recurrence given in Eq. (123). These expressions then relate 
the effective intercept parameter of the n-particle correlation function to the effective 
intercept parameter of the two-particle correlation function, and could be thus checked 
experimentally. For that type of test, a few warnings should be made: The intercept 
parameters A* jn (P) explicitly depend on the momentum of the particles in the region 
where all the n particles have approximately the same momentum. Thus averaging over 
the transverse mass and the rapidity of the particles to improve statistics is in principle 
not allowed. Further, the effective intercept parameter A #)Jl (P) differs from the exact 
analytical value of the n-particle correlation function which is n! in our picture, due to the 
fact that we have neglected possible partial coherence of the particle emitting source. Our 
motivation for ignoring possible partial coherence was to find a relatively simple limiting 
case, and also we are aware of the fact that the pure quantum statistical relationship 
between the second order and the higher order correlation functions has been shown not 



to be consistent with the available UAl data [12|]. In case of heavy ion reactions, resonance 
halo seems to be describing the drop of the effective intercept parameter from the value 



of 2 both in numerical simulations |22| and in analytical approaches [|Tj| 13fl . Finally, 
we stress that the values for the effective intercept parameters A*, n (P) may depend on 
the two-particle momentum resolution Q m i n , on the error distribution on the correlation 
function, as well as on the Gaussian or non-Gaussian structure for the Fourier-transformed 
Wigner-function of the core, S c (K,Ak). 

Note also that the general result for the correlation function in the core/halo model 
coincides with a particular limiting case of the expressions obtained by Cramer and Kadija 
for the correlation functions of the order 1, 2, 3, 4 and 5, namely, the case when particle 
mis-identification (contamination) is taken into account but the source is assumed to have 
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no partially coherent component. Thus our results can be applied also to the evaluation 
of the intercept parameter of higher order correlation functions for data contaminated by 
unidentified or errorneously identified particles. 

Prediction of higher order intercept parameters for NA44 data. Equations ( f24] - |2"8| ) can 
also be used to predict the intercept parameter of higher order correlation functions from 
the measured intercept parameter of the second order correlation function. The published 
NA44 correlation data indicate that A2,* — 0.52 ± 0.02 approximately independently of 
the transverse mass for both the low p t and the high p t data sample for pion pairs in a 
central S + Pb reaction at 200 A GeV. If we interpret this A2,* in the core/halo model, 
we obtain that f c = 0.72 is approximately independent of the transverse mass in the 
NA44 acceptance range. Thus, the core/halo model predicts the following measurable 
intercept parameters for the higher order correlation functions of identified pions in 200 
AGeV S + Pb reactions in the NA44 acceptance: 

A*, 3 (P) = 2.3 A*, 4 (P) = 8.6 A*, 5 (P) = 33.4 A*, 6 (P) = 148.0 (30) 

It should be emphasized that the NA44 data constitute a good sample to test the core/halo 
model on higher order correlation functions, given the experimentally observed momentum 



independence of the A 2 ,* parameter |3T 



In summary, we have evaluated the higher order correlation functions in the core/halo 
model analytically. A closed form of the correlation function of arbitrary high order is 
given by eq. fllTD in terms of the momentum-dependent core fractions and in terms of 
permutations that completely mix a set of indexes such that p. t ^ i. A recurrence relation 
has been found which allows for the evaluation of the measured intercept parameter of 
the n-particle correlation function for arbitrary high orders in a very efficient manner. 
We emphasized that the separation of the core from the halo is dependent on the relative 
momentum resolution of the experiment and pointed out a formal analogy between the 
correlation functions of the core/halo model and those of completely chaotic sources with 
particle mis-identification. Our results allow for a prediction of the higher order correlation 
functions if the basic building block, the amplitude s(i,j) is determined experimentally as 
a function of the mean and the relative momentum of the particle pair. If this quantity is 
real, it can be determined from a detailed analysis of the two-particle correlation function. 
If the imaginary part of s(i,j) is not negligible, a simultaneous analysis of second and 
third order correlation functions is necessary to extract the building block of higher order 
correlations. 
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